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ABSTRACT 

Most stars in galactic disks are believed to be born as a member of star clusters or 
associations. Star clusters formed in disks are disrupted due to the tidal stripping 
and the evolution of star clusters themselves, and as a results new stars are supplied 
to the galactic disks. We performed TV-body simulations of star clusters in galactic 
disks, in which both star clusters and galactic disks are modeled as TV-body ("live") 
systems, and as a consequence the disks form transient and recurrent spiral arms. In 
such no n- steady spiral arms, star clusters migrate radially due to the interaction with 
spiral arms. We found that the migration timescale is a few hundreds Myr and that 
the angular momentum changes of star clusters are at most ~ 50% in 1 Gyr. Radial 
migration of star clusters to the inner region of galaxies results in a fast disruption of 
the star clusters because of a stronger tidal field in the inner region of the galaxy. This 
effect is not negligible for the disruption timescale of star clusters in galactic disks. 
Stars stripped from clusters form tidal tails which spread over 1-2 kpc. While the 
spatial distribution of tidal tails change in a complicated way due to the non-steady 
spiral arms, the velocity distribution conserve well even if the tidal tails are located 
at a few kpc from their parent clusters. Tidal tails of clusters in galactic disks might 
be detected using velocity plots. 
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dynamics — methods: N-body simulations 



1 INTRODUCTION 

Star clusters are one of the fundamental building blocks 
of galact ic disks because m ost stars are formed in star 
clusters (|Lada fc Ladal l2003h . In disk galaxies, new clus- 
ters born in galactic disks travel in their host disks 
experiencing disruptions and supply new stars to the 
disks. The disruption of clusters is caused by the tidal 
force from their host galaxy and also the internal evolu- 
tion of star clusters themselves such as dynamical evolu- 
tion, ma ss loss due to the stellar evolution, and gas ex- 
pulsi o n (iFukushige &; Heggid 2000l; Baumgardt &; Makinol 



20031: iDehnen et al.1 12004 ICapuzzo Dolcetta et alJ l2005fc 
Lamers et alJ 120051 : iBaumgardt et al.1 120081 : iGieles etal] 
201 lh . 

Non-axisymmetric structures in gal actic disks, 
such as spiral arms (IGieles et al. I l2007h and bars 
([Berentzen Athanassoulal 2012f h have also been ex- 
pected to play important roles for the dynamical disruption 



of clusters. For example, IGieles et al.1 (|2007h investigated 
the effect of spiral-arm passages on the evolution of star 
clusters assumed to rotate in a fixed pattern speed, as 
in the stationary d ensity wave theory (|Lin fe Shul Il964l : 
iBertin fc Lin I Il996h . However, self-consistent simulations 
of galactic disks have shown that self-excited spiral arms 
are not stationary regardless of the existence of gas (or 
some kind of dis si patio n) (iSellwood &; Carlberg 1984 : 
I Sellwood fe Binnevl 1200% iBottemal 120031: iBaba et al I 



* E-mail: 

baba.j .ab@m.titech.ac.jp( JB) 



fujii@strw.leidenuniv.nl(MSF); 



120091: iFuiii et al.1 120111: ISellwood! liuTTl : IWada et al.1 120111 : 
iGrand et 811 2012: Ba ba et alJ l2012). Such non-steady spiral 
arms do not have a si ngle pattern speed but roughly follow 
the galactic rotation (jWada et al.ll201ll : [Grand et al.ll2012h . 
Therefore, these arms scatter st ars everywhere in the disk 
by the co-rotation resonance (pellwood fc Binne v 2002; 
IGrand et all 120121 : IBaba et al.1 boi2h . Non-steadv spiral 
arms change the gravitational fields around star clusters 
chaotically rather than periodically as in the stationary 
density waves. 

In order to know the dynamical evolution of star clus- 
ters in galactic disks with non- steady spiral arms, we need 
to model both star clusters and galactic disks as iV-body 
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("live") systems. Such self-consistent iV-body simulations 
are technically more difficult than those with rigid poten- 
tial disks because the dynamical timescale of star clusters is 
much shorter than that of galactic disks and a large num- 
ber of particles a re required for the modeling of the disks 
([Fuiii et al.ll201lh . We solve this probl em using a direct-tree 
hybrid code, Bridge ([Fuiii et al.ll2QQ7h . 

In this letter, we perform self-consistent iV-body simu- 
lations of star clusters in live disks using Bridge and demon- 
strate that the angular-momentum exchange between star 
clusters and spiral arms causes the radial migration of star 
clusters of a few kpc from their initial galacto-centric radii. 
The migration timescale is shorter than the galactic rota- 
tion timescale, i.e., a few hundred Myr. The radial migration 
causes the tidal disruption of star clusters bringing them to 
closer to the galactic center. Star clusters lose their mass 
in their perigalacticon passage, and their tidal tails spread 
over a few kpc. We also find that tidal-tail stars stay close 
to their parent clusters in their velocity space even if they 
are already a few kpc from the parent cluster. The tidal tail 
of young clusters might be detectable using their velocities. 
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2 iV-BODY SIMULATIONS 

We performed a series of iV-body simulations of star clus- 
ters embedded in a live galactic disk with spiral arms. We 
modeled both the disk and the clusters as iV-body sys- 
tems, but the halo is modeled as a potential. We s et up 
the disks following models used in iFuiii et al.l (|201lh . We 
adopted an exponential disk model with a total disk mass 
of 3.2 x 1O 1O M with 3 x 10 6 (3M) particles and as a con- 
sequence the mass of a disk-particle is ~ 1O 4 M0. The scale 
radius and scale height of the disk are 3.4 kpc and 0.34 
kpc, respectiv ely. For the dark mat ter halo, we adopted the 
NFW model (jNavarro et al.lll997T ) with the concentration 
parameter of the halo, c — 10. The virial radius and the 
mass of the halo are 122 kpc and 6.4 x 10 n M©. We mod- 
eled star clusters as a Kin g model wit h the dimensionless 
central potential Wo = 3 (|Kind Il966h . We adopted a to- 
tal cluster mass of 10 5 M© and a half-mass radius of 8 pc, 
and therefore the tidal radius is 30 pc. Our model is simi- 
lar t o young massive clusters in M51 and M82 (see Figure 
9 in IPortegies Zwart et al.1 (|2010h ) rather than those in in 
the Milky Way disk, which are more compact and therefore 
would be tidally disrupted less than our model. Included the 
stellar evolution, however, the Milky- Way clusters might ex- 
pand a factor of 5-10 i n the first 10 Myr as seen in obser- 
vations (see Figure 8 in IPortegies Zwart et al ] (|2010h ). We 
used 8192 (8k) equal-mass particles for the cluster. 

We first integrated only the disk up to 5 Gyr, in which 
self-excited spiral arms fully developed from the initial Pois- 
son noise due to the swing amplification. Then, we detected 
dense regions in the disk using the procedure below and put 
star clusters in the dense regions assuming that they are 
born there. We detected disk particles whose eighth-nearest- 
neighbor position is closer than 70% of the Jacobi radius of 
a cluster with 1O 5 M0. We chose the densest position from 
the detected locations and rejected other candidates within 
five Jacobi radii from the selected one in order to avoid that 
star clusters initially collide with each other. Repeating this 
procedure, we chose 97 positions in a spiral arm in 4-10 



Figure 1. Snapshots of the initial condition and at 0.25, 0.5, and 
1 Gyr. Gray scale shows the column density of the disk, and color 
points indicate cluster particles. Since these are only twelve colors 
but 97 clusters, we plot several clusters with the same color. The 
color points are as large as the cluster size to make the points 
visible. More than 80% of cluster-particles are still bound at 1 
Gyr. 

kpc from the galactic center. The initial positions of the 
star clusters are shown in the top left panel of Figure [T] We 
adopted the center-of-mass velocity of the eight neighbors 
as the cluster velocity. 

The initial radial, azimuthal, and vertical velocity dis- 
persion among the star clusters are 10.2, 7.7, and 7.5 km s _1 , 
respectively. Although these values are smaller than the 
mean of the disk stars (19, 14, and 14 km s _1 at 8 kpc 
respectively), star clusters would have rather smaller veloc- 
ity dispersion than those of old stars if we assume that star 
clusters form from giant molecular clouds as is observed in 
the Milky Way (the Geneva- Copenha gen Survey of the Solar 
Neighborhood; iHolmberg et al.ll2009l ). 

The simulatio ns are performed using a direct-tree hy- 
brid code, Bridge (|Fuiii et al.ll2007h . In Bridge, only the in- 
ner motion of star clusters are integrated using a direct N- 
body code and the other intera ctions are integrated using 
a tree code ([Barnes &; Hut| [T986V We a dopted a sixth-order 
Herm ite scheme for the direct method ([Nitadori fe Makinol 
2008) without any softening and with an accuracy parame- 
ter of 0.9. For the disk particles, the gravitational potential 
is softened using Plummer softening with a length of 10 pc. 
We adopted an opening angle of 0.4 with the center-of-mass 
approximation and a time step of 0.29 Myr for the tree code. 



3 RESULTS 

3.1 Radial migration of star clusters 

Figure [T] shows the time evolution of the star clusters and 
disk. Although each spiral arm looks like a single coher- 
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ent structure in each snapshot, it is transient and recurrent. 
They are wound up due to the differential rotation of the 
galactic disk and break up into multiple segments with a 
few kpc-scale, but the segments reconnect and form new 
long coherent arms (see online material movie). 

The orbits of star clusters in such transient spiral arms 
are complicated rather than a simple epicyclic motion. The 
motion of star clusters in spiral arms is similar to that o f 
gas and stars in the disk as shown in Wada etal ] (|201lh : 
I Grand et al ] (|2012h . Clusters tend to stay in or close to spi- 
ral arms and migrate along the arms. One of the orbits (the 
distance from the galactic center) obtained from our simula- 
tion is shown in the top left panel of Figure [2] (see also online 
material movie). The cluster is initially located at around 9 
kpc and migrates inward down to around 6.5 kpc losing its 
angular momentum, but it migrates outward up again to 8.5 
kpc (see the top left panel in Figure [2]) . 

Such radial migration is caused by the dynamical inter- 
action between star clusters and spiral arms. The bottom 
left panel of Figure [2] shows time evolution of the azimuthal 
force from the disk on the star cluster. The cluster loses 
its angular momentum, when it is moving ahead of a spiral 
arm (see top right panel of Figure [2j in which spiral arms 
are moving from right to left). On the other hand, the star 
cluster gains angular momentum, when it is moving behind 
a spiral arm (see middle right panel of Figure [2]). The an- 
gular momentum of clusters does not change when clusters 
are located at just the middle of two arms (see bottom right 
panel of figure [2]). Thus, star clusters in non-steady disks 
lose or gain angular momentum and as a result migrate a 
few kpc from their initial positions. 

Such angular-momentum changes occur for all star 
clusters. The clusters in our simulation lost or gained at 
most ~ 50% of their initial angular momenta within 1 
Gyr, which corresponds to radial migration of a few kpc. 
The angular-momentum change in the clusters is similar 
to that of stars investigated in simulations o f stellar disks 
(ISellwood Binnevll2QQ2l : iRoskar et al.ll2008l : iGrand et all 
l2012l : lBaba et alJ2012h . Since star clusters in transient spiral 
arms migrate a few kpc in their galactic rotation timescale, 
i.e., a few hundred Myr, open clusters in the Galactic disk 
older than ^100 Myr are expected to have already migrated 
from their initial galact centric radii. 

This rapid migration of star clusters may also change 
our understanding of the evolution of the Galactic disk. 
Indeed open clusters in the Milky Way disk have been 
used as tracers for t he dynamical and chemical evolution 
of the Galactic disk (jFriel fe Janedll993l : Ichen et al.ll2QQ3l : 
iMagrini et al.ll2009L e.g.,). In these studies, the spiral arms 
are considered to be stationary density waves, but these re- 
sults might change in non-steady arms. 

Another important effect of the radial migration is the 
tidal disruption of star clusters due to the smaller Jacobi 
radii at smaller distance from the galactic center. In our 
models it changes from 64 pc at 8 kpc to 45 pc at 4 kpc. 
In Figure [3] we see a clear correlation between the bound 
mass at the end of our simulation (1.5 Gyr) and the mini- 
mum perigalacticon distance of the clusters. In our simula- 
tion, clusters lose masses mainly during their perigalacticon 
passages, and therefore the mass loss becomes larger when 
they migrates inward. H owever, internal he ating due to the 
passage of spiral arms ([Gieles et al. 1 120071 ) does not work 
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Figure 3. The relation between the bound mass at 1.5 Gyr 
normalized by the initial mass and the minimum perigalacticon 
distance. 



efficiently in the case of transient arms, because both star 
clusters and spiral arms are always corotating. 



3.2 Tidal tails 

Stars tidally stripped from their parent clusters form tidal 
tails spreading over a few kpc (see Figure [T]). Their shapes 
are more complicated than those in simple spherical external 
potentials and change in time. When star clusters migrate 
from apogalacticon to perigalacticon, they move along a spi- 
ral arm, and their tidal tails also elongate along the spiral 
arm. When clusters are in the other phases (e.g., perigalac- 
ticon or apogalacticon), the tails come closer to the parent 
cluster again even if they are no longer bound to the clus- 
ter (see top right panel of Figure [2]). This behavior is sim- 
ilar to that in the case in axis ymmetric external potentials 
(ICapuzzo Dolcetta et aDl2005h . 

Figure U shows spacial (left) and velocity (right) distri- 
butions of two clusters. The cluster shown in the top panels 
is at its perigalacticon passage, and the one in the bottom 
panels is moving from its perigalacticon to apogalacticon. 
The cluster shown in the top panels is the same as that 
shown in Figure [2] In spacial distribution plots (left pan- 
els in Figure 2]), cluster particles (crosses) show tidal tails 
spreading over a few kpc which are close to the orbits (black 
doted curves) or the circular orbits at the positions of the 
clusters (black dashed curve). In spite of such a large spa- 
tial distribution, we find that the tidal-tail stars still remain 
close to their parent clusters in velocity space. In the right 
panels, we plot the position of cluster stars in vr-v^ space as 
crosses. Colors in both panels indicate the velocity deviation 
from the parent clusters; blue, red, green, and black indicate 
<2, 2-5, 5-7, and >7 km s _1 , respectively. We also plot disk 
particles within 1 kpc from the cluster center (cyan points). 
In contrast to the cluster particles, the disk particles have 
a larger distribution in the velocity spaces. In spatial plots 
we plot disk particle which are within 1 kpc from the clus- 
ter center and have velocity deviation of <5 km s _1 in the 
velocity space (cyan points). If we detect star cluster tails 
using only their velocities, these stars would be detected as 
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Figure 2. Left column: Evolution of a star cluster. Top: radial distance from the galactic center. Middle: angular momentum. Bottom: 
the force in azimuthal direction from the disk on the cluster normalized by the force in radial direction. Right column: position of the 
cluster particles (red points) and the center of mass (blue cross) and the column density of spiral arms in the R — <fi plane at the times 
indicated by dashed lines in the left column. From top to bottom the star cluster is located inter, forward, and behind a spiral arm, 
respectively. 



contamination. They are located randomly in spatial plots, 
while star cluster particles distribute in tidal tails. 



4 SUMMARY 

We performed a series of iV-body simulations of star clusters 
in live stellar disks with multiple spiral arms. In these sim- 
ulations, both galactic disks and star clusters are modeled 
as TV-body systems and integrated self-consistently. Our re- 
sults show that star clusters migrate radially a few kpc in the 
time scale of their orbital period in the disk (a few 100 Myr) 
because of the angular-momentum exchange with transient 
spiral arms. The angular-momentum change of the clusters 
is at most ~ 50% of the initial angular momentum within 1 
Gyr. 

In the case of non- steady transient spiral arms, the ra- 
dial migration strongly affects the tidal disruption of star 
clusters because star clusters lose more mass when they ap- 



proach the galactic center due to the smaller Jacobi radii. 
This kind of disruption mechanism does not appear in sta- 
tionary density waves. The heating due to the non-steady 
(corotating) spiral-arm passage would not be as strong as 
that by the density- wave spiral arms, because the adiabatic 
change of the energy due to the slow passages of spiral arms 
suppresses the heating per passage, and the number of spi- 
ral p assages is qu it e few (see the corotation case in Figure 
8 in Gi eles et al. I (|2007h ). Furthermore, the radial migra- 
tion of star clusters can carry stars far from their original 
orbital radii and finally the distribution of stars would be 
much wider than those in t he case of the density waves 
(iBland-Hawthorn et ai]l20ld ). 

With transient spiral arms, star clusters and their tidal 
tails tend to stay in or close to spiral arms. The shape of 
tidal tails of clusters change in a complicated way in time 
compared to those in a smooth tidal field like a halo poten- 
tial. When a star cluster approach the galactic center, star 
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Figure 4. Distribution of cluster particles and disk particles 
around clusters in x-y plane (left) and vr-v^ phase (right). The 
cluster shown in the top panels is the same one shown in Figure [2] 
at a time of 1.75 Gyr, when the cluster is located at its perigalac- 
ticon. The bottom panels show a cluster, which is moving from its 
perigalacticon to apogalacticon at a time of 1 Gyr. Crosses and 
points indicate cluster and disk particles, respectively. In the left 
panels, colors indicate the velocity shown in right panels. Black 
dotted and dashed curves show orbits obtained from the simu- 
lation and circular orbits at the distance of the cluster from the 
galactic center, respectively. Yellow circles show the 1 and 2 kpc 
radius from the center-of-mass position of the cluster. In the right 
panels, blue, red, green, and black crosses indicate cluster parti- 
cles with <2, 2-5, 5-7, and >7 km s _1 from the mean velocity 
of the cluster particles. Cyan points indicate disk particles which 
are located within 1 kpc from the cluster center. Disk particles 
which are located within 1 kpc and have velocities <2 and 2-5 
km s — 1 from the mean velocity of the cluster are shown in the 
left panels as blue and red points, respectively. 

clusters move along a spiral arm as is the case of stars and 
gas moving in spiral arms. In this phase, the tidal tails also 
spread along the spiral arm. During the apogalacticon pas- 
sage, on the other hand, the tidal tails are compressed and 
distribute around 1 kpc from the cluster even though they 
are unbound. The tidal tails of clusters might be detectable 
even if they spread over a few kpc, because the tidal-tail 
stars still remain very close to the cluster in velocity space 
after they become unbound. If we know the velocity of stars, 
we might be able to detect the tidal tails of star clusters in 
the Galactic disk using a future astrometry such as Gaia 
and JASMINE. 
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